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ABSTRACT 

We use the AdS/CFT correspondence to argue that large dyonic black holes in anti-de Sitter 
spacetime are dual to stationary solutions of the equations of relativistic magnetohydrodynamics on 
the conformal boundary of AdS. The dyonic Kerr-Newman-AdS4 solution corresponds to a charged 
diamagnetic fluid not subject to any net Lorentz force, due to orthogonal magnetic and electric 
fields compensating each other. The conserved charges, stress tensor and R-current of the fluid are 
shown to be in exact agreement with the corresponding quantities of the black hole. Furthermore, 
we obtain stationary solutions of the Navier-Stokes equations in four dimensions, which yield pre- 
dictions for (yet to be constructed) charged rotating black strings in AdSs carrying nonvanishing 
momentum along the string. Finally, we consider Scherk-Schwarz reduced AdS gravity on a circle. 
In this theory, large black holes and black strings are dual to lumps of deconfined plasma of the 
associated CFT. We analyze the effects that a magnetic field introduces in the Rayleigh-Plateau 
instability of a plasma tube, which is holographically dual to the Gregory-Laflamme instability of 
a magnetically charged black string. 
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1 Introduction 

The AdS/CFT correspondence (cf. [1] for a review) has provided us with a powerful tool to get 
insight into the strong coupling dynamics of certain field theories by studying classical gravity. 
In the long wavelength limit, it is reasonable to expect that these interacting field theories admit 
an effective description in terms of hydrodynamics. One can thus conversely use the equations of 
fluid mechanics (which, in simple contexts, are much easier to solve than the full set of Einstein's 
equations) in order to make predictions for the gravity side. This was done in [2], where the fluid 
configurations dual to (yet to be discovered) black rings in (Scherk-Schwarz compactified) AdS were 
constructed^ . 

More recently, it was shown in [4] (cf. also [5] for the four-dimensional case) that the equations 
of hydrodynamics (i. e., the Navier-Stokes equations) can also be derived directly from Einstein's 
equations with negative cosmological constant, without making any use of string theory, from which 
these ideas originally emerged. 

The correspondence between gravity in asymptotically AdS spaces and fluid mechanics has 
interesting consequences: For instance, it is well known that under certain conditions fluids are 
affected by the so-called Rayleigh-Plateau instability (responsible e. g. for the pinch-off of thin 
water jets from kitchen taps), and one might ask what the gravity dual of (the relativistic analogue 

^For a discussion of plasma lumps dual to black saturns in AdS cf. [3]. 
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of) this instability is. This point was studied recently in [6], where it was argued that the gravity 
dual is the Gregory-Lafiamme instability [7] . 

In this paper, we shall consider a more general setting in which the charged fluid moves in 
external magnetic fields. The fluid is thus described by the eqiiations of magnetohydro dynamics 
(MHD). Such a generalization is of interest for several reasons: First, there are two possible fall-off 
conditions for a bulk U(l) gauge field: A normalizable mode that corresponds to a VEV of the dual 
operator (an R-current), and a non-normalizable one corresponding to the application of an exter- 
nal gauge field, that can be thought of as an electromagnetic field. (2+l)-dimensional field theories 
deformed in this way, that are dual to magnetically charged AdS4 black holes, have become fashion- 
able recently in the context of a possible holographic realization of condensed matter phenomena 
like superconductivity [8,9], the Hall effect [10] and the Nernst effect [11]. In the long wavelength 
regime a quantum field theory in presence of external gauge fields is expected to be described by 
MHD. Second, solving the MHD equations (which, under certain symmetry assumptions, might be 
much easier than solving the full set of Einstein-Maxwell equations) can be helpful for constructing 
new magnetically charged black hole solutions in AdS. Furthermore, the phase structure of such 
AdS black holes can be studied in a simplified setting using magnetohydrodynamics. 

We should stress that the magnetic field entering the MHD equations is nondynamical, so that 
there are no Maxwell equations on the boundary. Note however that, in the same way in which 
one can promote the CFT metric to a dynamical field [12], it should be possible to add dynamics 
also to the magnetic field, although we will not try to do this here. 

The remainder of this paper is organized as follows: In section 2, we review general aspects of 
MHD, like the equation of state and the stress tensor of the fluid under the influence of external 
magnetic fields. Moreover, we show how the MHD equations emerge from Einstein-Maxwell- AdS 
gravity. In the following section, we generalize the results of [13] to the case of nonvanishing 
magnetic fields. We first consider static dyonic black holes in AdS4. In an appropriate limit 
(when the horizon radius r/j is much larger than the AdS curvature radius ^), these black holes are 
effectively described by magnetohydrodynamics. Conformal invariance and extensivity imply that 
the grandcanonical partition function has to take the form 

^lnZg, = T^hiC/T,B/T^), (1.1) 

where C denotes a chemical potential conjugate to the U(l) electric charge of the fluid, B is the 
applied magnetic field, and V and T represent the volume and the overall temperature of the fiuid 
respectively. We then show that, using the function h that can be read off from the partition 
function of static black holes as an input into the MHD equations, we can exactly reproduce the 
conserved charges, boundary stress tensor and R-currents associated to the rotating dyonic Kerr- 
Newman-AdS solutions. Their thermodynamics is very simple: It is summarized by the partition 
function 

with H = 1 — 6<j^£^, and lo is the angular velocity of the black hole. 

After that, we solve the Navier-Stokes equations on M x S*^ x S^, IR x S*^ x or R x S"^ x H'^. 
This yields predictions for (yet to be constructed) charged rotating black strings in AdSs with 
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momentum along the string. 

In section 4 we consider perturbations around plasma tube solutions of the MHD equations 
on R X X S^, that are dual to magnetic black strings on Scherk-Schwarz compactified AdSe- 
These plasma tubes suffer from the long wavelength Raylcigh-Platcau instability, which is shown 
to be weakened by the presence of a magnetic background. This is exactly what one expects from 
the known studies of the Gregory-Laflamme instability for magnetically charged black strings in 
asymptotically flat space, to which the considered SS-compactified AdS black strings are similar. 

Throughout this paper we use calligraphic letters T, V, i3, 5, . . . to indicate local thermodynamic 
quantities, whereas T,V,B, S, . . . refer to the whole fluid configuration, n and ( are local and global 
chemical potentials respectively. 

Note added: While this paper was in preparation, ref. [14] appeared that partially overlaps with 
our section 2.3. 

2 Magnetohydrodynamics 
2.1 Equation of state 

We would like to know how the equation of state V = p/{d—l) of a conformal fluid in d dimensions 
changes in presence of a magnetic field B. To this end, consider the grandcanonical potential 

i^ = £-rS- ijlR = ^{T,V,ijl,B). (2.1) 

It follows from conformal invariance and extensivity that 

$ = -VT'^h{u, h) , (2.2) 

where we defined v = jjl/T and b = B/T'^. Note that /x and B have mass dimension one and two 
respectively, so that ly and b are dimensionless. Equation (2.2) defines the function h{u,b). Prom 
(2.2) it is also clear that 

$(Ar, A^-'^V, A/x, X^B) = A$(r, V, n, B) . (2.3) 

If we derive this with respect to A, set A = 1 and use 

df~~ ' 5V~ ' 'd^ 
where M is the magnetization density, we obtain 

-sr - (1 - d)rv - 7^/x - 2vmb = $ . (2.5) 

Using the definition (2.1) of this yields the equation of state 

p={d-l)V - 2BM . (2.6) 

We will see below that large dyonic black holes in AdS4 indeed have a grandcanonical potential of 
the form (2.2), with d = 3. 



= -7^, 



dB 



= -VM, 



(2.4) 
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2.2 Stress tensor of the fluid 



We want to determine now the general form of the stress tensor for our fluid, and we will do this 
in a derivative expansion in the fluid velocity u^, following [15]. The zero order term corresponds 
to a perfect fluid (under the influence of external electromagnetic fields F^^, labelled by the index 
/), with stress tensor, charge currents and entropy current given respectively by (cf. e. g. [11]) 

T^''' = pu^'u" + Vn'"' -MfF^'^x, J'; + VaM''/ = riu^' , 4 = su^" , (2.7) 

where we have introduced the projection tensor 

l^ni^ = g^^-^ + uf^u" (2.8) 

on the directions orthogonal to u^. p = p{T, p^ ,B'^) is the rest frame energy density, rj = 
ri{T, p/\B^) and s = s(T, p^ ,B'^) arc the rest frame charge- and entropy densities, while p^ 
and B"^ denote the chemical potentials and magnetic fields respectively. The polarization tensor 
M.j'^ is defined by 

Mf = ■ (2.9) 

Note that the microscopic currents (polarization currents) are given by 

whereas Jj represents the total current density, so that the combination Jj + Va-M'j'^ appearing in 
(2.7) is the macroscopic or transport current. Notice also that the extra contribution —Ai^^F^^\ 
to the stress tensor in (2.7) comes from the coupling of the polarization to the electromagnetic field. 

The first subleading order in the derivative expansion of the MHD equations describes dissipative 
phenomena like viscosity and diffusion. Lorentz invariance and the requirement that the entropy is 
non-decreasing determine the form of the stress tensor and the currents. Let us introduce diffusion 
currents g^, 

+ VaMT = riu^ + , (2.11) 

with the constraint ti^Qj = 0, meaning that the diffusive process is purely spatial according to an 
observer comoving with the fluid clement. Next wc introduce the heat flux and the viscous stress 
tensor, decomposed in a symmetric traceless part t^^ and a trace r, such that 

Ti^v = + ^)n'^i^ + pu^u"" - MfF^"^ + q^'u'' + q^'u^ + t^"" . (2.12) 

Again these quantities are purely spatial, and verify the constraints 

Ui,q^ = Q, ii;.T'^" = 0. (2.13) 

Note that the new fields , g'^, r and r^^ are first order in Vn. Finally the entropy flux is a linear 
combination of all the available vectors that are at most first order in Vu. Therefore 

= su^ + ^q^ - X^q^j , (2.14) 
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where the scalars /3 and are functions of the thermodynamic parameters. 

To determine the most general form allowed for the newly introduced quantities we impose the 
second law of thermodynamics. To have a non-decreasing total entropy the four-divergence of the 
entropy flux must be positive, 

V^J|>0. (2.15) 
Projecting the MHD equations on we obtain 

u^V.T'^^ = n^jrF^, (2.16) 

which can be cast, with some straightforward manipulations, in the form^ 

-lu^^MfV^Fl^ - u^Fl,{r: + V^MT) , (2.17) 

with the expansion = V^n'^. 
We have the Euler relation 

p + r = sT + n^ri, (2.18) 
which implies the Gibbs-Duhem relation 

dV = sdT + ndu^ + MidB^ , (2.19) 

with the magnetization densities Mi- This yields 

sV^r = V^:P - r/V^M^ - MiVf^B' . (2.20) 

The entropy density is a function of p, rj and B"^ , so that its gradient reads 

(2.21) 

Then 

^.r, = i.„.v,p - t^u^v^r. + ^u^v,B^ + ( ^ + 7 ) 

+13^ fit + 'f^^ll - A' V,.,;' - qlV^X' . (2.22) 
The conservation of the charge currents Jj = yields the relation^ 

u^'V^ri = -V^q'^-rii^, (2.23) 

^To get this, one has to use Mj^ = XijF'^^^ ^ with the susceptibihties xu ~ XJit as well as the Bianchi identities 
for F'. 

^Note that the polarization current is separately conserved, V^VaA^'^'* = 0. 
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that, together with (2.17), can be used to put the divergence of the entropy flux in the form 

-q^ (V^A^ - ^F^X^ - - ^T^'^V.u, . (2.24) 

Note that, in order to get (2.24), we assumed that the grandcanonical potential <I> depends on 
only through the invariant F^^^F^f^'^ (no summation over /), and we defined (B^)'^ to be given (up 
to a prefactor) by F^^F^'^''. This leads to 

-^MfV^^Fl^ = -MiV^B' . (2.25) 

We ensure that the right hand side of (2.24) is positive or vanishing by requiring that each term 
be positive or vanishing. This can be easily achieved by choosing 

^=1, ^' = Y- (2.26) 



Then, the diffusion currents are given by 



(2.27) 



where we have introduced the projector 11^'^ to ensure that the diffusion currents be spacelike, and 
the diffusion matrix Djj is positive definite. Moreover, we set 

T = -CV^u'^ , (2.28) 

where C ^ is the bulk viscosity, and the heat flux is given by 

with K> the thermal conductivity of the fluid. Finally, 

r'^'^ = -na"" , (2.30) 
with 7] > the shear viscosity of the fluid and a^'^ the shear tensor, 

= 1 (U'^'VaU^' + W^VaU") - ^m'^l? . (2.31) 

2 3 

With these definitions, we find 

Vm^5^ = ^ + + i^-') + ^ ^ ' (2.32) 

which is positive by construction. Notice finally that the stress tensor in (2.7) is traceless on account 
of the equation of state (2.6). 
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2.3 MHD equations from gravity 

In this section, we consider Einstein-Maxwell- AdS gravity in D dimensions and show that at leading 
order the boundary theory dual to a charged Z)-dimensional AdS black hole reduces to magnetohy- 
drodynamics md = D—1 dimensions in the long wavelength sector. We follow the usual procedure 
to foliate asymptotically the spacetime with timelike hypersurfaces and regularize the action 
by adding appropriate boundary counterterms. Then, in the r ^ oo limit, the Brown- York stress 
tensor has a finite limit, the renormalized holographic stress tensor of the dual CFT. The pro- 
jection of the Einstein- Maxwell equations on S,. then shows that this stress tensor together with 
the boundary gauge fields satisfy the equations of magnetohydrodynamics on the boundary. We 
can interpret this as the leading contribution in the derivative expansion to the boundary MHD 
equations describing the long wavelength sector of gravity. 

The Einstein- Maxwell action with negative cosmological constant A = —{D — 2){D — l)/2^^ 
reads 



/ = J <Px^g {R-2A- J'mnJ''''') + ^cs + ^ d'^xV^ K + I,t, (2. 



33) 



where Ics is the Chern-Simons term, present in the D = 5 case. M, N, . . . are bulk indices, whereas 
^,1/, . . . refer to the boundary. In this subsection (and only here), we refer to the bulk gauge field 
as J^MN to distinguish it from the boundary gauge field. In all other sections of the paper we do 
not need such a distinction. denotes the boundary hypersurface r = R with outward pointing 
unit normal , and induced metric 

hMN = QMN - riMfiN ■ (2.34) 

K is the trace of the extrinsic curvature defined by (with Vm the bulk covariant derivative) 

Kmn = hM^VpUN ■ (2.35) 

Finally, Id are the usual boundary counterterms needed to obtain a finite action in the limit 
i? — > oo. We shall not need the precise form of let, but only the fact that their variation with 
respect to the metric is divergence-free, 

V'^-^=0. (2.36) 

Here we have defined as the induced covariant derivative on S/j. Also, notice that additional 
terms can be added to let to handle the logarithmic divergences appearing for odd Z), corresponding 
to the Weyl anomaly of the dual CFT. We can use the Fefferman-Graham expansion to write the 
metric of any asymptotically AdS spacetime near spatial infinity in the form 



ds^ = —dr^ + r^5^,.(r, x)dx^'dx'', (2.37) 



where 



9.^ = + + • • • + + '^ + (^) . (2.38) 
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The coefficients gf^} and hi>li}, depend only on the boundary coordinates and the coefficient , 
related to the holographic Weyl anomaly, is present for odd D only. The normal vector to the 
constant r hypersurfaces is n^dA = r/£dr and the induced metric is h^i, = r'^g^i/\^^^- Then, the 
conformal boundary metric of AdS d is obtained by taking the limit 

7^^ = lim hij,y . (2.39) 

R—*oo jx 

Let us decompose the gauge field in the orthogonal component and its projection Fmn on S^, 

J^MN = Fmn + ^{JMnN - JNriM) ■, (2.40) 

where 

jM ^ , Fmn = hM^'hN^J'pQ ■ (2.41) 

The analysis of the asymptotic behavior of vector gauge fields [16, 17] shows that Ff^i, has a finite 
R ^ CO limit, while the current J*^ goes to zero like r^'^. Therefore, the renormalized background 
gauge field in the dual field theory, and the i?-symmetry current read respectively 

F^^= lim F^^, J^'= lim . (2.42) 



Now, by varying / with respect to the bulk metric and gauge field, we obtain the equations of 
motion^ 

Fmn = Gmn + ^Qmn — Tmn = , T^^^mn = , (2.43) 
where the stress tensor of the gauge field is given by 

Tmn = "^^M^^PN — ■^^'^Qmn ■ (2.44) 

Then, from the projection 

EMNU^h^P = (2.45) 
of the Einstein equations, we obtain, using Gauss-Codazzi 

Vm {K^n - S^nK) = FnmJ^ . (2.46) 

We recognize in the term in parenthesis the Brown- York boundary stress tensor, which diverges as 
we take the oo limit, but will give a finite limit - the holographic stress tensor T'^'^ - once we 
take into account the counterterm contribution [18], 



-7 l^'T^p = lim h^^^ - K,K + ^ , (2.47) 



SIct \ 

Since the counterterm contribution is divergence-free by (2.36), we can add its divergence to the 
left hand side of (2.46) and, after multiplying the equation by \/— 7i, obtain a finite i? — oo limit 
that reads 

V^T'^. = F.^J'' , (2.48) 



*In D — 5 one has to take into account that there is a Chern-Simons term and the Maxwell equations read 
d*J^ + -^T AJ^ = 0. 
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where is the covariant derivative induced on the conformal boundary with metric On the 
other hand, using the definition (2.41) it follows that V^J'^ = 0, that becomes, in the i? — oo 
limit, the conservation law for the i?-current, 

V^J^ = 0. (2.49) 

Note that equations (2.48) and (2.49) are simply the Ward identities associated to the bulk diffeo- 
morphism and gauge invariance respectively [19]. 

These conservation equations for the boundary stress-energy tensor and i?-current become the 
equations for magnetohydrodynamics in the background field F^i, when a large black hole is present 
in the bulk. Indeed, as noted in [20], stationary Kerr-AdS black holes have a holographic stress 
tensor that assumes the perfect fiuid form. We will show in the next section that this is still the 
case for dyonic Kerr-Newman-AdS4 black holes, and we believe it to be true in any dimension. 
This dual stress tensor, of the perfect fluid form, yields, when combined with equations (2.48) and 
(2.49), the equations of MHD. Then, if one perturbs these stationary solutions with long wavelength 
disturbances, in the spirit of [4], the horizon can still be decomposed into patches that tubewise 
extend to the boundary which are approximated by boosted pieces of black branes. Therefore, 
to leading order in the derivative expansion, large magnetically charged black holes are dual to a 
magnetohydrodynamic theory. We shall check this explicitely in the next section for dyonic AdS4 
black holes. Higher orders in the perturbation theory produce higher-derivative dissipative terms 
in the stress tensor. A detailed analysis of this long wavelength sector of gravity, with a complete 
proof of the duality with MHD and the computation of the dual stress tensor up to third order in 
the derivative expansion has been performed by Hansen and Kraus in the AdS4 case and appeared 
during the last stages of preparation of this manuscript [14] . 



3 AdS black holes and black strings from MHD 

3.1 Static dyonic black holes in AdS4 

The equations of motion following from the Einstein-Maxwell action with negative cosmological 
constant A = — 3£~^, 

admit the static dyonic black hole solutions 



ds 



2 - -V{r)dt^ + + r\de^ + siefdcp') 



where 



and 



V{r) 



V{r) 



+ k 



2m ql + 



+ 



n2 ' 



s{e) = { 



sin 6 , k = I , 
1, k = 0, 
sinh 6 , k = —1 . 



(3.2) 
(3.3) 
(3.4) 
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The horizon is thus S'^ {k = 1), {k = 0) or (k = —1). m, qe and qm denote the mass 
parameter, electric and magnetic charge respectively. The one-form gauge potential reads 



r 



= qm j s{e)de . (3.5) 



The strength of the magnetic field B in the dual CFT can be obtained (up to rescaling by powers 
of t) by taking r — >■ oo in the expression for the bulk field strength. This leads to B = qmj^ ■ The 
electric charge density Pe of the state in the field theory is given by = (J*), where is the 
R-current that can be computed as follows. On-shell we have for the variation of the action with 
respect to the gauge potential 
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6A^ 



Miv = y d^xdM{V^F''''6AN) = / d^xV^nMF^''6AN , (3.6) 



where h denotes the induced metric on the boundary, and n is the outward pointing unit normal 
to the boundary. One has thus in the limit of large r 

^ = J-^ (37) 
SAt 47rG r3 ' ^ ^ ^ 

To get the CFT R-current, one has to rescale this by r^/l'^, so that 



where we used the AdS/CFT dictionary 

_J__V2N^ 
levrG ~ 247r£2 



(3.9) 



Note that the result (3.8) was obtained for k = m [10]. In order for the potential to be regular 
at the horizon r = rh, At must vanish there. This requires that wc add the pure gauge term 
[qg/rfii)dt to A/i. This term is non-normalizable and has the dual interpretation of adding a 
chemical potential for the electric charge, ( = qe/vh^, to the field theory [10]. The R-charge is 
R = PeV, with the spatial volume^ 

V = e'^ J ded(ps{e) . (3.10) 

This yields 

The entropy S and energy E of the black hole read 



A,_V2N^ 



order to get a finite volume in the cases k = and k = —1, one has to compactify the horizon to a torus or a 
higher genus Riemann surface respectively. 
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E = 



m 



In what follows, we shall be interested in a magnetohydrodynamical description of the above black 
holes. Like in [13], we may estimate the mean free path for the fluid as Zmfp ~ v/p^ where rj is the 
shear viscosity and p is the energy density of the fluid. For fluids described by a gravitational dual, 
one has r] = s/Att, where s is the entropy density [21]. Consequently, Imfp ~ s/inp, and an MHD 
description will be valid if this value is much smaller than the radius of the S"^ or the curvature 
radius^, s/Airp <^ £. Using the expressions for S and E, this implies 



(3.14) 



A sufficient condition for (3.14) to hold is r^/i 3> 1, i. e. , for large black holes. In this case we can 
neglect the ^-contribution to the energy, and the thermodynamic fundamental relation becomes 



E(SVRB)-^^^^(^^y^' 



2 ^2jv3y2 ^ 
^ 1852 



One easily checks that 



where 



dE_ ^ _^ 



T = 



1 

in 



^2 



(3.15) 



(3.16) 



(3.17) 



is the Hawking temperature of the black hole. From (3.15), wc get the grandcanonical potential 



^T,V,C,B) = E-TS-CR 
with the function h given by 



h 



247r 



-VT^h{C/T, B/T"^ 



_3_ ( B_ 
H VT2 



where H{(/T, B /T'^) = rh/{l'^T) is determined by the fourth order equation 

3iJ^ - AttH^ - {C/TfH'^ - {B/T'^f = , 

that follows from (3.17). Note that $ has indeed the form (2.2), as it must be. 

We now want to obtain the dyonic AdS4 black holes from MHD on ]R x ^2 , . 
The metric on the conformal boundary is given by 

ds'^ = -dt^ + £2(d02 + S{efd(f)^) , 

so that the only nonzero Christoffel symbols are 

r<t> -r<t> - ^'^^^ 



or 



= -s{e)s'{e) , 



s{e) 



(3.18) 



(3.19) 



(3.20) 



xiJ2. 



(3.21) 



(3.22) 



In the (compactified) fc = case, Zmfp must be much smaller than the length of any torus cycle. 
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For stationary, translationally invariant and axisymmetric configurations one has dtT^^ = d(j)T^'^ = 
(we also assume dtF^"" = d^Pi^'' = 0), and the MHD equations V^T'^'" = F'' ^J^' become thus 

a^T^* + ^T^* = F^Ji" , (3.23) 

deT^'^ + ^r^^ - SS'T^'^ = F^Jf" , (3.24) 

doT^'l' + ££_r^^ = i?^^^ . (3.25) 
The macroscopic electric charge current and the entropy current are given respectively by 

JiLcr. = J'' + '^aM'"' = ru^, 4 = su^, (3.26) 

where = 'y{l,v) is the 3-velocity of the fluid, r denotes the electric charge density and s is the 
rest frame entropy density. Both currents are conserved, 

V;.J^acr.=V;.4 = 0. (3.27) 

As there are no dissipative terms in the charge- and entropy currents, we have for the entropy 

s = Jcfx^/^4 = j ded(i)i^s{e)s-i , (3.28) 

and for the electric charge 

R = j (fx^Ji,^,, = J d9d<t>eS{e)r^. (3.29) 

The Kilhng vectors of interest are dt (energy E) and (angular momentum J on the S"^, or 
H"^). The conserved charge related to a Kilhng vector k is proportional to / d^Xy/^T^ ^k^, and 
hence 



E = J ded(t)fS{e)T*\ (3.30) 

J = j d9d(t){^S^{9)T*'l' . (3.31) 

The shear tensor cr^^, heat flux q'^ and diffusion current (7^ must vanish on any stationary solution 
of fluid dynamics. The requirement a'^'^ = means that the fluid motion should be just a rigid 
rotation. By an SO (3) transformation^ we can can choose this rotation such that the 3-velocity of 
the fluid is u*^ = (u*, u^, u^) = 7(1, 0, a;) for some constant uj. Prom 'u'^n^ = —1 and 7 = (1 — ■u^)~^/^ 
one obtains then v"^ = u'^ l'^ S'^ (9) . 

The equilibrium fluid flow is symmetric under translations of t and (p, so that all thermodynamic 
quantities depend only on 9. We now evaluate the expansion, acceleration, shear tensor, heat flux 
(2.29) and diffusion current (2.27), with the result 

^ = 0, a" = {0,-S{9)S'{9)jW,0), ai"' = 0, (3.32) 



^Strictly speaking, an applied electromagnetic field breaks S0(3) invariance, so that our choice of u'^ implies 
consistency conditions on F^i, (cf. (3.35) below). 
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,0 



(3.33) 



(3.34) 



where k and D denote the thermal conductivity and the diffusion coefficient respectively, T is the 
local temperature and /j, the local chemical potential. The requirement that and vanish 
implies that 



r 



T7: 



H4> 



0. 



7 



0, 



(3.35) 



with r constant. If Fg„u" = 0, the last equ. is solved by /x = Ti/, where is constant. The 
conditions (3.35) determine all the thermodynamic quantities as a function of the coordinate 9. We 
now want to shew that this configuration solves the equations of magnetohydrodynamics. To this 
end, we first notice that the dissipative part of the energy-momentum tensor vanishes once (3.35) 
is imposed, so that all nonzero contributions to the stress tensor result from the perfect fiuid part 
plus interaction with the external electromagnetic field, and read 

f p + v'^V {p + V)u \ 

\ (p + 7')a; pu? + Vi-'^S-'^' 
where V denotes the local pressure and 7]^^ is given by 



(3.36) 



rp\XV 



( M'^Fte 


-M^'f'Fi 







te 



-l-^M^'^Fte + 




-M^^Ft 


-e-^S-^{9)M^'l'Fe^ J 



(3.37) 



Eqns. (3.23) and (3.25) imply F^^J^^ = F'f'^J^^ = 0, which are automatically satisfied if the sus- 
ceptibility X does not depend on t and (f), which we assume in the following. The only nontrivial 
equation of motion (3.24) becomes 



(3.38) 



Using the Gibbs-Duhem relation (2.19) as well as (2.25), equ. (3.38) can be cast into the form 



7S 



d_ 



+ r7 



d 



rFe^ui^ 



(3.39) 



which is automatically solved using (3.35). Therefore, the rigidly rotating configurations are sta- 
tionary solutions of the magnctohydrodynamic equations. Moreover, as the diffusion current and 
the heat flux vanish, the only nonzero contributions to the transport R-charge- and entropy currents 
come from the perfect fluid pieces (3.26), which are easily seen to be conserved as well. 
Prom the grandcanonical potential (2.2) we find 



T^h{u,b), M=T 



dh 



3h 



dh 
/ — 

dv 



2b 



dh 
db 



p = 2{V - MB) = 2T^ (h - b^'j 



r = T- 



.dh 
dv 



(3.40) 
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In the nonrotating case a; = 0, we have 7 = 1, Fq^u'^ = Fq^u'I' = Fq^^ui = 0, B = B = const., 
and thus (3.35) imphes that the local temperature T as well as u and b are constant. We can then 
easily compute the energy, entropy, R-charge and magnetization, with the result 



db J ' du ^ 

, _V(3.-.|^-..|). M...|. (3.41) 

It is straightforward to verify^ that the Hawking temperature of the black hole is given by T = r, 
and the chemical potential associated to the R-charge is C = tz^- Using the relations 

dh V2Ar3/2 ^Ar3/2 ^ 

ai; = ^^^^' db = -^^H^ ^'-''^ 

following from (3.19), (3.20), one easily shows that E, S and R coincide with the corresponding 
expressions (3.15), (3.12) and (3.11) for the static black hole. 

In the next subsection we will show that, using the function h determined from the static AdS4 
black hole as an input into the MHD equations, one can exactly reproduce the thermodynamics, 
boundary stress tensor and R-current of the rotating dyonic Kerr-Newman-AdS4 solution. 

3.2 Kerr-Newman-AdS4 

The Kerr-Newman-AdS4 solution to the equations of motion following from (3.1) is given by^ 

ds^ = -^[dt - ^d^f + ^dr^ + + ^^[adt - ^#]^ , (3.43) 

A = -S^ldt - + '-^[adt - , (3.44) 

where ^ ^ 

Ar = ir^ + a^)(^l + ^^^-2mr + ql + ql, ^ = 1-^, (3.45) 

= r^ + cos^ 6* , = 1 - cos^ Q , (3.46) 

and a is a rotation parameter. 

The thermodynamical phase structure and stability of these black holes was studied in detail 
in [22]. In order to get the boundary geometry da^, one takes r = const. — 00, and then rescales 
ds^ by i^jr^. This leads to 

,0 asin^e*,,,, £2 jj^2^ 

da^ = -\dt ^ — d^f + —de^ + ^ d(j)^ , (3.47) 

which is a rotating Einstein universe. In the limit of large r, the U(l) field strength behaves as 

F ^ ^sm9d0 Adcp. (3.48) 



One can express dE — rdS — rudR + MdB in terms of dr. du, dh and check that it vanishes. 
^For the time being, we restrict to the case of spherical horizons. The generahzation to fe = 0, — 1 is straightforward. 
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As before, the expectation value of the CFT R-current can be computed by varying the action 
on-shell with respect to the gauge potential A^^ and subsequently rescaling by . This yields 

posin Q ^ ^ ,„ 

= J^ = ' •^^ = 0- (3-49) 

The chemical potential C, of the dual field theory is given by [22] 

C-^ = ■^'vX Ir— »oo ~ -^ZX \r=rh = 9 ' (3.50) 

where X = dt + ^ndcfi denotes the null generator of the horizon and A^, is the electric part of the 
vector potential. The angular velocity of the horizon reads 

= (3.51) 

In what follows, we shall be interested in the limit of large black holes, Vh S> when a hydrody- 
namical description is valid. As a? < i'^ (otherwise Ag can become negative), we have then also 
S> o^, and the chemical potential (3.50) reduces to 

The electric charge of the black hole (the field theory R-charge) is obtained by computing the flux 
of the electromagnetic field strength at infinity, with the result^^ 

R=^. (3.53) 
The temperature, Bekenstein-Hawking entropy, energy and angular momentum are given by [22] 
rh f, a\ 3rl + ql + 1 /Sr^ ql + q^ 



r 



h 



s = ^ ' '\ (3.55) 



-Kj rj + g^) _ W2iV3/2r2 



E = 



m x/2Ar3/2 |-^3 ^ ^2 ^ ^2 1 



J = aE, (3.57) 

where in the last steps we took the large black hole limit. Let us write the metric (3.43) in the 
ADM form 

ds'' = -N^dt^ + ^dr^ + ^d9^ + ^^|^(# - ndtf , (3.58) 

Ar Ag p^^^ 

where 

T? = {r'^ + a'^fAe-a^ArS\n^e, (3.59) 



""Note that we rescaled Rhy & factor of £ and C by £ ^ with respect to [22]. 
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and the lapse function A'^ and the angular velocity are defined by 

^ = ^[Ae(r2 + a2)-A,] . (3.60) 
It was shown in [22] that the angular velocity u entering the thermodynamics is 

u; = nH-n^ = ^^^4^^, (3.61) 

which boils down to 

LU « a/f (3.62) 

for large black holes. 

By applying the implicit coordinate transformation [23]^^ 

T = t, ^ = (f)+^t, ycos@ = rcos9, = ^[r^A^ + sin^ 6*] (3.63) 

to the solution (3.43), one gets a different slicing such that the boundary geometry is that of a 
static Einstein universe, 

da^ = -dT^ + £2(de2 + sin2 ed$2) . (3.64) 

Note that (3.64) is related to (3.47) by a diffeomorphism plus a Weyl rescaling: On the boundary 
r — GO (which implies y oo), the coordinate transformation (3.63) reduces to 

Qj sm 

T = t, ^ = (f>+—t, sine = ^=. (3.65) 



Applying this to (3.64) and subsequently rescaling with e^'^ = Ag/H yields (3.47). The two bound- 
ary metrics are thus conformally related. For the five-dimensional Kerr-AdS black hole, this was 
first noticed in [24]. 

In the new coordinates, the Faraday tensor on the boundary becomes 

/ 2 \ -3/2 

F ^ qmU-^sin^ej sin 9 dG A - -^dr) . (3.66) 

Due to the nonvanishing component Fq^i we have now also an external electric field. This is not 
surprising; it results from the boost in (3.65). 

In order to proceed, we need the thermodynamic variable B. We already stated in section 2.2 
that is proportional to F^^F^^ . To fix the prefactor, note that for the field strength (3.66), one 
gets 



^^We apologize for using the same symbol $ for an angular coordinate and the grandcanonical potential, but the 
meaning should be clear from the context. 
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which reduces to = 2g^/^*^ = 2i^B'^ in the static case a = 0. We have thus^^ 



\/F2 

B = . (3.68) 

The grandcanonical potential $(T, V, /Lt,H) is given by (3.18), with h{i',b) defined in (3.19), (3.20), 
and 

6 = — = ^ ^ . 3.69 

One can then compute the polarization tensor M^^, with the result 

1 

M''" = -y,^=xF''\ (3.70) 

with the susceptibility 

_ I dh_ V2N^/^ 

Note that x is always negative, so that our fluid is diamagnetic. Moreover, for large (local) tem- 
peratures T (keeping fi and B finite), (3.20) gives H{iJ,/T,B/T^) —>■ 47r/3, and thus 

i. e., the susceptibility obeys a diamagnetic Curie law with Curie constant C. 

Let us now consider the MHD equations V ^T^'^ = ^J^^ with the field strength (3.66). One 
easily shows that F'^ ^u^^ = for z/ = T, <5, and that F'^f^u^ is proportional to a; — a/i'^. Let us 
assume that this vanishes, i. e. , that the angular velocity of the fluid is given by a; = a/£^, as is 
suggested by (3.62). Then the transport current ru^^ is orthogonal to the electromagnetic field, and 
(3.35) implies that u = ji/T is constant. Moreover, from (3.69) and (3.67) we have 

which is constant by virtue of (3.35). Using this, it is straightforward to shew that the polarization 
current J^ipj. = — Vo-TW^^ is proportional to in this case, and hence orthogonal to F^ ^ as 
well. This means that F^ ^^J^^ = 0; due to orthogonal magnetic and electric fields compensating 
each other, there is no net Lorentz force acting on the fiuid, and the MHD equations boil down 
to VjuT'*^ = 0. It is interesting to consider this situation from the point of view of a reference 
frame which is moving with the charged fluid (the rotating Einstein universe (3.47)). In this 
reference frame the fluid is static and therefore not subject to the magnetic force. Since the net 
electromagnetic force, which in this frame consists only of the electric force, must be zero, the 
electric field must vanish in the moving reference frame. This is indeed the case, as can be seen 
from (3.48). 



^^Up to the prefactor i, that stems from the fact that the boundary Faraday tensor is obtained from the corre- 
sponding bulk quantity by projecting on the boundary and then rescaling with , in fiat spacetime this would just 
be the statement that = 2{B^ — £^), and we defined to be what is usually called — £^ . 



17 



It would be very interesting to see to which kind of black holes the fluid configurations with uj 
different from a/i'^ correspond to. They might have to do with a nontrivial nut parameter, but we 
shall leave this point for future investigations. 

One can now proceed to compute the conserved charges. Prom (3.28) - (3.31) one obtains 

47r£V2 / . dh ^,dh\ „ AttI'^t'^ dh 



S = — R-—^^, (3.74) 
E = ^[h-i-\, , = , (3.75) 



together with the magnetization 

M 



The temperature T and chemical potential ^ are given by T = r and = tu respectively, while 
the intensive variable conjugate to M is 5 = r'^b/E. To show this, one can express dE — rdS — 
TudR + MdB — udJ in terms of dr, du, db, dw and check that it vanishes. Note that T, ( and 
B are distinct from T, /j, and B. While the former quantities are asssociated to the whole fluid 
configuration, the latter are local thermodynamic properties of the fluid that vary on the boundary 
manifold. For the grandcanonical partition function one finds 

1 VT'^h(C/T.BE/T^) 
\nZg, = -±.{E-TS-u;J-CR)= ^ , (3.77) 

which means that InZgc of the rotating fluid is obtained simply by dividing the grandcanonical 
partition function of the static fluid by the universal factor S = 1 — u^i^, and replacing B by BE. 
This generalizes the results of [13] to nonvanishing magnetic fields. 

We must now compare the fluid charges with the corresponding black hole quantities. To do 
this, we first note that for large KNAdS black holes, the temperature has the same dependence on 
r/i as in the static case (compare (3.54) with (3.17)), so that the function H = rfi/{i'^T) is again 
determined by (3.20). Using this together with (3.19) and (3.42), one easily shows that the fluid 
charges R, S, E, J exactly coincide with the expressions for the black hole given in (3.53)-(3.57). 
This coincidence extends also to the energy-momentum tensor and the R-current: The holographic 
stress tensor of the KNAdS black hole was determined in [22], and reads 

„ 2m ^ 2ma sin^ 9 / „ x 



Th = ^^^^2 ' Tee - , (3.79) 

and all other components vanishing. This corresponds to the boundary geometry (3.47), i. e. , 
to the rotating Einstein universe. To get the stress tensor for the static Einstein universe (3.64), 
recall that the latter is related to the former by a diffeomorphism plus a Weyl rescaling, da^ = 
e^'^da^, with e^'^ = Ag/E. Under a Weyl rescaling, the energy-momentum tensor transforms as 
T^^ = e-('^-2)'^r^^ (see e. g. [13]), which yields (taking d = 3) 
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= 7 (37 - 2) , Tee = , (3.81) 

which is easily shown to exactly coincide with the fluid stress tensor (3.36). 

The R-current (3.49) transforms as = e~'^^J^^ under a Weyl rescahng [13]. This gives for 
the R-current corresponding to the static boundary geometry 

This is again equal to the transport current Jmacr. = ru^ of the fluid. 
3.3 Black strings in AdSs 

We now solve the Navier-Stokes equations on M x S"^ x 5^, M x 5^ x or M x 5^ x i^^. This 
will yield predictions for black strings in AdSs, for which at present only partial results are known 
[25-31]. Since the function h{i',b) entering (2.2) is unknown for d = 4, 6 7^ 0, we shall consider 
hydrodynamics rather than MHD, i. e. , we will restrict to the case of zero electromagnetic fields 
in this section. 

The metric entering the MHD equations is 

ds^ = -df + dz^ + l^ide"^ + S{e fd(l?) , (3.83) 

where S{9) was given in (3.4). The only nonzero Christoffel symbols arc the ones in (3.22). For 
stationary, translationally invariant and axisymmetric configurations one has dtT^^^ = dzT^^ = 
d(f)T^'^ = 0, and the Navier-Stokes equations reduce to (3.23)-(3.25) (with F^i, = 0) plus the 
additional condition ^ 

deT^' + ^T^' = . (3.84) 

The entropy and the R-charges read 

S = Jd^x^/^Jl = J dzded(l)fS{9)s-f, (3.85) 

Ri = Jd^x^/^Jj = J dzd9d(t)e'^S{9)rn. (3.86) 

Now the Killing vectors of interest are dt (energy E), d^p (angular momentum J on the S"^, or 
H^), and dz (momentum p along the string). The associated conserved charges are 

E = J dzd9d(l)e^S{9)T*^ , (3.87) 

J = J dzd9d(t>i^S^i9)T^'^ , (3.88) 

p = J dzd9d(l)£'^S{9)T*^ . (3.89) 

The shear tensor a^'^ , heat flux and diffusion currents qj must vanish on any stationary solution 
of fluid dynamics. The requirement a^'^ = means that the fluid motion should be just a rigid 
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rotation. By an SO (3) transformation we can choose this rotation such that the 4- velocity of the 
fluid is = {u*,u^,u^,u'^) = 7(1, wi, 0, 0^2) for some constants (jJi,(jJ2- Prom n'^u^ = —1 and 
7 = (1 - ?;2)-i/2 one obtains then v'^ = + £'^u^S'^{9). 

The equilibrium fluid flow is symmetric under translations of t, z and (f), so that all thermo- 
dynamic quantities depend only on 9. We now evaluate the expansion, acceleration, shear tensor, 
heat flux and diffusion current, with the result 



q^ = -r'Kj{0,0,i-(-],0] , (3.91) 



Q'^ = -r'Du{0,0,^(^^),0) . (3.92) 



^ = 0, = {0,(},-S{9)S'{9)-f'^ulO), a'"' = 0, (3.90) 

d9 Vt, 
^ //r 
d9\r 

The requirement that and q'^ vanish implies that 

r = r7 , / = Tu^ , (3.93) 

with r and constant. The conditions (3.93) determine all the thermodynamic quantities as a 
function of the coordinate 9. We now want to shew that this configuration solves the Navier-Stokes 
equations. To this end, we proceed as in [13], and first notice that the dissipative part of the 
energy-momentum tensor vanishes once (3.93) is imposed, so that all nonzero contributions to the 
stress tensor result from the perfect fluid part, and read 

/ p + v'^V {p + V)lji {p + V)u)2 

{p + V)ui pul + V{l-g-'^ulS'^{9)) {p + V)ojiU2 

7"V^ 

\{p + r)uj2 {p + r)uxuj2 pul + rg'^s-^{9){i-ujl) J 



1' 



The only nontrivial equation of motion (3.24) becomes 



dV ^ d In 7 , , X 

Using the Gibbs-Duhem relation (2.19), (3.94) can be cast into the form 

which is automatically solved using (3.93). Therefore, the rigidly rotating configurations are sta- 
tionary solutions of the equations of fluid dynamics. 
Prom the grandcanonical potential 

$ = -VT^h{u^) (3.96) 

one obtains 

p = W = 2,T*h{u^), n = T^hi{iy^), s = T^{Ah-v^hi), (3.97) 
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with 

dh 

hi = -Pf. (3.98) 

We saw in the previous section that for large AdS4 black holes, the function h is insensitive to the 
curvature parameter k. The same is to be expected for large black strings in AdSs, so that we can 
infer the function h{u^) from the A; = case, which is simply the AdSs black hole with fiat horizon. 
This yields [13] 

where the auxiliary parameters are related to the by 

The conserved charges corresponding to our fluid configurations are 

E = Lfr^hiu') (4Ki,6K, u;2) " ^^1,4^, c^s)) , (3.101) 

p = ALfT'^uih{u^)Ki^6{ui, a;2) , J = 4L^V^a;2/i(i^^)K3,6(a;i, 0^2) , (3.102) 
S = Lfr^ {Ah - v^hi) ^1,4(^1, a;2) , Ri = LfT^hi{u'^)Ki^4{ooi,uj2) , (3.103) 
where L is the length of the parametrized by z, and we have defined the integrals 

Km,n{^i,u;2) = j S^{9)^^d(t>d9 . (3.104) 

These integrals, which can be performed in terms of elementary functions, satisfy the relations 

^"f"""'" = nUlKm,n+2 , ^"5'"'" = nfu2Km+2,n+2 , (3.105) 

OUl OU)2 

and 

-£'^UJ2Km,n = Km-2,n-2 " (1 " ^l)-f^m-2,n • (3.106) 

The chemical potentials corresponding to the rotating/boosted fluid solutions are defined by 



S,J,p,Rk 



One easily verifies that^^ 

T = T, W = CUl , Q = CU2 , = T/J,^ . (3.107) 

Using this, we find for the grandcanonical partition function 

In Zgc = -^ {E -TS -wp-nj -C^Rj) = Lt^T^h{u^)Ki^i{uJi,ui2) ■ (3.108) 



^^As before, one can express dE — rdS — Tv^dRi — LOidp — uj2dJ in terms of dr, dv^ , dwi, dw2 and check that it 
vanishes. In order to show this, one has to use the relations (3.105) and (3.106). 
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In the k = case, the integration is trivial and Zgc reads 

~ {l-ul-l-'ulf ' 



InZ, 



(3.109) 



where V is the volume of the three-torus. 
For A; = 1, we have 



arctan 



+ 



(3.110) 



where V = AttM'^ is the volume of the x S^. Note that this reduces, for loi = 0, to 



In Zgc = ^VT^h{u^) 



arctan 



+ 



eL02 



(3.111) 



while for u;2 = it becomes 



\lVZgc=-VT^h{v') 



+ 1 



(3.112) 



We see that the formula correctly reproduces the 1/(1 — u;^) behaviour characteristic of the rotation 
on spheres and the 1/(1 — w^)^ behaviour of the boost along a circle (compare also with the results 
of [13] for Kerr-AdS black holes and with [23]). 

Our results of this section yield predictions for rotating black strings in AdSs with momentum 
along the string, that carry three electric charges (more precisely, there are electric charge densities 
along the string). While the general solutions of this type are still to be discovered, there are some 
partial numerical results available, namely for static charged black strings with all three charges 
equal [29]. It would be very interesting to compare the thermodynamics of these solutions with 
our hydrodynamical predictions. Unfortunately, the gravity solutions have been constructed only 
numerically, so that the comparison involves a fair amount of numerical work that goes beyond the 
scope of this paper. 



4 Gregory-Laflamme on magnetized plasma tubes 

4.1 Description of the problem. Scherk-Schwarz compact ificat ion of a CFT 

It was observed in [2,32,33] that, in the long wavelength regime, d-dimensional fluid dynamics is also 
an effective theory describing the Scherk-Schwarz (SS) compactification of a (d + l)-dimensional 
CFT. This theory is dual to SS compactification of AdSrf-|-2 gravity. At finite temperature the 
theory has two vacuum solutions, namely i) the black branc solution that describes the dcconfined 
phase of the dual gauge theory and ii) the AdS soliton that describes the confined phase [32]. In the 
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neighborhood of the critical phase transition temperature the two phases can cohabit in equihbrium 
with a domain wall interface. Black holes in SS compactified AdS(i+2 are then dual to deconfined 
plasma configurations immersed in the confined phase. The phase diagram of plasma equilibrium 
solutions includes plasma balls, pinched plasma balls, plasma rings and plasma tubes [2,6,33]. 
These plasma solutions are dual to SS compactified AdS black objects whose explicit solution is 
not yet known. However, quite interestingly, these phase diagrams reveal to have properties that 
are quite similar (although with some differences linked to the AdS nature of the system) to the 
known phase diagram solutions of asymptotically flat black holes. In [6] we have argued that this is 
not a coincidence since asymptotically fiat black holes may indeed admit a fluid description in the 
limit of large number of dimensions. It is therefore important to explore the properties of plasma 
objects since they might provide a lighthouse to understand better even asymptotically flat black 
holes. We will return to this discussion in the end of this section. 

In a previous paper [6], we have shown that plasma tubes suffer from the long wavelength 
Rayleigh-Plateau instability that makes them pinch-off when their length is bigger than their ra- 
dius. We have further observed that these solutions are the natural duals of SS AdS black strings 
and that the plasma instability is dual to the Gregory-Laflamme instability of a black string. The 
aforementioned analysis was constructed on the top of a series of studies [34-39] where the remark- 
able analogy between fluid tubes and black strings was observed and analyzed. After [2,33], [6] 
promoted this curious analogy to an actual duality. 

In this section we want to go a step forward and discuss the effect that a magnetic background 
introduces in the Rayleigh-Plateau instability of a plasma tube. This will provide solid predictions 
for the features of the dual magnetic black strings, associated Gregory-Laflamme instability and 
phase diagram of associated magnetically charged solutions. 

We will take our boundary geometry to be x x 5^ , described by the metric 

In this background, we consider fluid configurations that arc uniform plasma tubes with topology 
X D^, i.e., a disk extended along the periodic z-direction and immersed in the vacuum or 
confined phase. We also take them to be translationally invariant along the SS direction (not 
represented in (4.1)). Furthermore, we assume that there is a constant magnetic field along the z 
direction. The dynamics of this plasma configuration is governed by the equations of relativistic 
magnetohydrodynamics (MHD) in (3 + l)-dimensions. As mentioned above, in the long wavelength 
regime, MHD is an effective theory describing the SS compactification of a (4 + l)-dimensional 
CFT. Moreover, in the dual gravity description these plasma tube configurations correspond to 
black strings in a background that asymptotes to a SS compactification of AdSe. The horizon 
topology of these strings is x S^. The reason being that on the interface between the plasma 
phase and the vacuum, the SS circle shrinks to zero size as we enter in the holographic direction 
and reach the horizon radius. 
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4.2 MHD for Scherk-Schwarz plasmas 

In this subsection we find the MHD equations that describe the long wavelength magnetohydro- 
dynamic limit of a Scherk-Schwarz compactification of a (d + l)-dimensional CFT. In the end we 
specialize to the d = 4 case. 

Neglecting subleading dissipation and diffusion contributions (in particular, in equilibrium these 
contributions vanish) the energy-momentum tensor of the fluid in the presence of an external 
electromagnetic field is the sum of the perfect fluid, Maxwell interaction and boundary surface 
contributions,^^ 

T;:,i = {P + 'P)u^u'' + rg^'\ 
T.I^ = -M^^F\, 

TZy = -<^hndf\6{f). (4.2) 

Here, u^^ is the fluid velocity, p, V and a are the density, pressure and surface tension of the fluid, 
F^i, is the external Maxwell field and M.^^ is the associated polarization tensor. The fluid boundary 
is defined by f{x^) = 0, it has unit normal = d^f /\df\, and h^'^ = g^" — n^n^ is the projector 
onto the boundary. 

The MHD equations follow from the covariant divergence of the stress tensor, 

V^T^"" = J,,F^''' , (4.3) 

where Ji, is the charge current defined in (2.7). Prom these equations we can derive the relativistic 

continuity, Navier-Stokes and Young-Laplace equations. In the absence of an electromagnetic field 
this derivation was presented in detail in [6]. With the extra Maxwell contribution the derivation 
follows similarly and we get (using footnote 2) 

ui'd^p + {p + r)V^u'' + ^u'^M'^^V^F^p = , (4.4) 

(p + r)u''V,u'' = -U^"' (v,r - ^M'^f^V.F^p^ + ru.Fi''' , (4.5) 

n< - n> = aiC , n = r- M^'^F^c^n^n" , K = hi^^V^n" , (4.6) 

where K is the boundary's extrinsic curvature and H is the generalized fluid pressure. It is the 
sum of the usual pressure V and the magnetic pressure (that is due to the Lorentz force exerted by 
the electromagnetic field on the polarization current). n< — n> is the generalized pressure jump 
when we cross the boundary from the exterior (with pressure n>) into the interior (with pressure 
n<). In the derivation of (4.6), the constraint that the fluid velocity must be orthogonal to the 
boundary normal is used (this guarantees that the fluid is confined inside the boundary) , 

uf'n^ = . (4.7) 



'^''The perfect fluid and Maxwefl interaction terms were already introduced in (2.7). In section 2, where we studied 
the MHD equations of a CFT on M'*, these were the only non-dissipative contributions. Here we want to study 

them in the case where wc have a SS compactification of a (4 + l)-dimcnsional CFT. Then there is also a surface 
contribution. It arises because in this case we have an interface between the plasma and the vacuum phases. 
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Equations (4.4)-(4.7) constitute the set of relativistic MHD equations governing the fluid dy- 
namics. A particular solution of these equations describes a static uniform plasma tube in the 
background (4.1) with radius Rq and extended along the z-direction, and with a constant magnetic 
field H(o) in the same direction. This solution is characterized by 



w 



(0)--^*' %)=^(0)K, ^ = ^(0), P = Pio)- (4.8) 



So, the non- vanishing Maxwell tensor components are Fr^ = —F^^ = B^^-^^J—g and F^j^F^*^ = 
2;B^pj. Wc use the subscript (0) to emphasize that this is a static unperturbed solution of MHD. 
In these conditions some terms in (4.4)-(4.6) simplify considerably. In particular, Ai'^I^S/ yFap = 
2X(o)V^;B(o) = 0, ru^Ff"" = and H = :P - M (0)8^0) (where M(q) = x^{o))- ^(0) is then a 
constant by the Navier-Stokes equations (4.5). The pressure ■P(o)) energy density /9(o) ^"^^ magnetic 
field ;B(o) are not independent. They arc related through the appropriate equation of state for a SS 
compactification of a (4 + l)-dimensional CFT. We derive this equation in the next subsection. It 
is given by (4.16) with d = 4. This equation of state together with the continuity equation (4.4) 
demands that p(o) is also a constant. 

In subsection 4.4 we will perturb equations (4.4)-(4.7) and study the Rayleigh-Plateau instability 
of a relativistic plasma tube in the presence of a magnetic field. 

4.3 Equation of state for a Scherk-Schwarz plasma 

In this subsection we discuss the grandcanonical potential and equation of state of a Scherk-Schwarz 
compactification of a (d + l)-dimensional CFT. The d = 4 result will be needed in the next 
subsection. 

Start with a. {d+ l)-dimensional CFT on M'^ x S^, where is the SS circle, in the presence of 
a conserved R-charge density r (and associated chemical potential fi) and of a constant magnetic 
field B. Conformal invariance and extensivity demand that the grandcanonical potential, 

^{T,V,fi,B) =£ -TS - fiTZ, (4.9) 

is given by 

$ = -Vr'^+^%,6), (4.10) 

where we defined u = fi/T and b = B/T"^. So, a CFT on M'^ x has the same thermodynamical 
potential as a CFT on M''+^ (that we analyzed in section 2.1-*^^). The compact SS direction does 
however introduce a Casimir or vacuum energy density po^^- Therefore, the equation of state of 
the CFT on M'* x is 

p = dV-2BM + po. (4.11) 

It follows from $ = ^{T,V, ij,,B) and (4.10) that the pressure of the conformal fiuid on M'^''"-'^ is 
given by 

= T'^+^h{iy,b). (4.12) 



V = 

dV 



'^^Note that in section 2.1 wc worked in d-dimcnsions while in this section we work in d-|- 1 dimensions. 



^In what follows, we shall assume that po is independent of v and b. 
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Replacing this in the equation of state (4.11) we obtain for the local temperature 

(p-po + 2BM\^^ (4.13) 

We can now do a dimensional reduction along the SS direction. This yields a field theory (that 
is not conformal) and which has a grandcanonical potential given by 

$,edCr,v,//,^) = £-rs-nn 

= v(^po-T'^+^h{u,b)) . (4.14) 

Upon this reduction, the local temperature of the plasma does not change and is still given by 
(4.13). Prom $red = *^red(^) V,/x,;S) One has, 

s = = '^'^ [{d + l)h{u, b) - ud,h{u, b) - 2bdbh{u, b)] , 



T''+'h{u,b)-po, 



dV 



1 5$red 



V dp 

M = ~^=r^-'d,h{u,b). (4.15) 
Replacing (4.13) in the expression (4.15) for the pressure we get the equation of state, 

r = ^ - + ^^P' + ^^-^ , (4.16) 

which is valid in or out of equilibrium. A similar replacement could be done for the other equations 
in (4.15). 

It is important to observe that the hydrodynamical equilibrium of the magnetized uniform 
tube also implies its thermodynamic equilibrium [6]. Indeed, the static condition d^V = that 
follows from the Navier-Stokes equation, together with the Gibbs-Duhem relation (2.19), written 
as VjyP = s Vj^T + r Vjy/x + Viy6, implies that the local temperature T is constant. Our uniform 
tube is static and thus the plasma temperature T equals the local temperature T. Demanding then 
that (4.13) is constant we find the relation (for constant c) 

p- po = cdh{iy,b) -2BM, (4.17) 

for an equilibrium static solution. The two last relations will be used later in the d = 4 case. 



4.4 Rayleigh-Plateau instability in a magnetic background 

We now address the stability of a uniform plasma tube with a constant magnetic field along the 
tube direction when we perturb it. The dynamics of the perturbations are dictated by the MHD 
equations subject to appropriate boundary conditions. Once perturbations settle in, the static 
plasma tube is taken away from thermal equilibrium and viscosity and diffusion effects start being 
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also active. Therefore, the energy-momentum tensor of the fluid includes now not only the perfect 
fluid, the Maxwell interaction and the boundary surface tension terms (4.2), but also a dissipative 
contribution. The uniform plasma tube is afflicted by the Rayleigh-Plateau instability [6]. Surface 
tension is the mechanism responsible for this instability. Viscosity and difi^usion play no role on the 
activation of the instability. Therefore, in our analysis we will neglect the dissipation contribution 
to the fluid stress tensor, and comment at the end on the effects it introduces. Our main aim is to 
find the effect that a magnetic field introduces on the instability, extending the analysis of [6]. 

Take a static uniform plasma tube with radius i?o in a constant magnetic field B^q^ along the 
tube z-direction. This unperturbed solution is described by (4.8) and by the equation of state 
(4.16), with cZ = 4. 

Consider now perturbations on this plasma tube. A generic perturbation is described as 

u^ = ufo)+<5n^, V = V^o)+SV, p = p(o)+5p, = ^fo) ' (4-18) 

where the perturbation on a quantity Q is generically denoted by 6Q. Note that the magnetic field 
is kept fixed. As there are no Maxwell equations on the boundary, it is treated as a non-dynamical 
parameter. Some of these perturbations are not independent. They are related by the equation of 
state (4.16), valid also out of equilibrium. Its perturbation yields the relation between the density 
and pressure perturbations, 

5p ~ ASV . (4.19) 

In deriving this, we made the assumption that S{BA4) = B^6x vanishes. This seems reasonable for 
small magnetic fields, or in a regime where the susceptibility varies only slowly as a function of the 
temperature and the other thermodynamic variables. We shall comment later on possible effects 
that a relaxation of this assumption might have. 

The perturbed state must satisfy the continuity and Navier-Stokes equations, (4.4) and (4.5). 
The Young-Laplace equation (4.6) and the constraint (4.7) provide boundary conditions for the 
perturbed problem. After eliminating the O*'^ order terms using the unperturbed MHD equations, 
the continuity and Navier-Stokes eqns. yield, up to first order in the perturbation, 

OrSu^ + ^Su^ + d^Su'f' + dt5u^ + dtSp = , (4.20) 

-n- P(0) + no) 

(P(o) + ^(0)) dtSuf" = -5WP - g'^'dJV + tSurF''^ + rSu^F'"'' = . (4.21) 

Since any perturbation can be written as a Fourier series, in the most general case we will 
consider a perturbation that disturbs the boundary surface of the plasma tube according to 

R = R{t, z, (j)) , R{t, z,(t)) = Ro + € e'^^e'''^+'""f' , e <^ Ro , (4.22) 

Positive CO describes an instability with wavenumber k. The possible unstable mode is axisymmetric 
if m = 0. Naturally, we look for perturbations of the fluid quantities that have the same form as 
the boundary disturbance, 

SQ{t, R, z, 4)) = SQ{R)e'^*e'''^+'"'^ , SQ = {Su", SV, 6p} , (4.23) 

which determines d^SQ for p = t, z, <p. 
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Solving (4.21) with perturbations (4.23) and (4.19) we find that the velocity perturbations at 
leading order are 



Su\R) = 
du^{R) = - 

Su'l'iR) = - 
5u^{R) = -i 



a^(P(0)+^(0)) (1 + ^ 
(P{o)+no))(l + § 



d5V{R) ^ imkA 



dR 



kR 



-1 



kA ddVjR) im 



k 



SViR) , 



w(P(o) +^(0)) 

where we defined the Alfven wavenumber of the system, 



kA = 



(0)' 



P(0) + 'P(o) 



(4.24) 



(4.25) 



Note that (4.24) satisfies the perturbed version of m^m^ = — 1, namely u^^^ 5u^ = 
Replacing (4.24) in the continuity equation (4.20) and using (4.19) we get 



dF5V{R) 1 dSVjR) 
dR"^ ^ R dR 



m 
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sr{R) = , 



r/^ = A;^ 1 + 1 + -f 



A;2 



A-2 



k^ 



(4.26) 



This is a modified Bessel equation. Its solutions are the modified Bessel functions of the first kind, 
ImivR), and second kind K^{r]R). Kjn{r)R) diverges as R~"^ as R ^ and we discard it (for 
m = 0, this solution would give a regular SV but Suf^ would still diverge; so we drop it also in this 
case). Therefore, the regular solution of (4.26) at the origin is 



dr{R) = AlmiriR) 



(4.27) 



where y4 is a constant to be fixed. Replacing the pressure perturbation in (4.24) yields for the 
radial component of the velocity, 



Su^{R) = - 



Ar] 



a;(p(o)+no))(l + &) 



(4.28) 



where Il{y) = dyl„{y). 

The solutions just found must satisfy a total of two appropriate boundary conditions. The 
first demands normal stress balance on the boundary. This means that the generalized pressure 
perturbation must also satisfy the third perturbed hydrodynamic equation, namely the one that 
follows from perturbing the Young-Laplace ewquation (4.6), 



BC I : 



K 



K\ 



Ro 



(4.29) 



Here, 5n|jjjj.y ~ 5'P|i,dry order e^. The reason being that the external magnetic field is 
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non-dynamical and thus fixed and the perturbation on the boundary normal is of order e^.^'' On 
the rhs we evaluate the expression at the perturbed boundary R = R{t, 4>, z) defined in (4.22) 
and we subtract the unperturbed contribution evaluated at i? = i?o- The extrinsic curvature is 
obtained from its definition, K = h^Vi^nf^, with the unit normal of the boundary (4.22), = 

\n\~^ (^—R^6^ + — R'^6^ — R'^5^. This condition fixes the ratio A/e to leading order as 

A a 



RllrnivRo) 



{k'^Rl + Lo'^Rl - 1 + m^) . (4.30) 



The second boundary condition is a kinematic condition requiring that the normal component of 
the fiuid velocity on the boundary satisfies the perturbed version of (4.7), u^^^^ + Su'^n^^^ = 0, 

where 6n^ = n^lj^^f. z (p) ~ '^m'^ unperturbed normal is njf'* = n^\j^^ = 6j}. This ensures that 

the velocity perturbation leaves the fiuid confined inside the boundary. This boundary condition 
then reads 

BC II : ^^i-^lbdry - ^ e e'^*e^"^^+'*=^ , (4.31) 

where the Ihs follows from (4.28) evaluated at the boundary. This boundary condition together 
with (4.30) yields the desired dispersion relation uj{k) for the Rayleigh-Plateau instability in the 
presence of a constant magnetic field along the plasma tube direction. 



a;2 = 



kRo Ii„.+i{i]Ro) , f-, , -kA 
+ m 1 + t — 



l + |f ) ' {1-m^-u^Rl-k'Rl) 



2X-I 



ImivRo) \ V 

(4.32) 

where we used the relation Il{y) = Iu+i{y) + jjluiv) and r] is defined in (4.26). p(o) + ^(o) is 
obtained from the equation of state (4.16) and (4.17) for d = 4. 

We conclude that the Rayleigh-Plateau instability is active (a; > 0) for wavenumbers k satisfying 
the condition 

(4.33) 

Since m is an integer this means that only axisymmetric modes (m = 0) are unstable. In a previous 
paper [6] we have shown that the condition (4.33) for the threshold wavenumber can be found by 
noting that unstable modes are those that reduce the potential energy of the system for fixed 
volume (equivalently, they are those that maximize the entropy for fixed conserved charges) . In the 
present system we have a competition between the surface tension and magnetic potential energies. 
However, to leading order 0{e) in the perturbation, the magnetic energy does not change. Indeed, 
the magnetic field is kept fixed without being dynamically perturbed, see (4.18), and the volume 
change that also contributes to the magnetic energy variation is subleading, i.e., 0{e^) (see footnote 
17), when compared with the 0{e) change in the surface tension energy. Therefore, only the surface 



^'^This point deserves a closer look. One has (5n|bdry = SV\bdTy — S (M^^F" (x'Tlfx'^u ) |bdry) S66 (4.6). Assum- 
ing 5x ~ 0, consistent with the approximation taken in (4.19), the magnetic contribution can be written as 



M^'"F"J{n^,n^) = M"" F\{nf^ Sn, + <5n^n[,">), where Sn^ = u^Irm.z) - 4°' and 4°^ = n^lfl^ = We then 
hsiveSiM^^F^^n^n^) = 2X(o)B(o)5nr- Now, using (4.22) we find that Srir = [(k'^ -u!^)/2 + m^ /{2Bi)]e'^ +0{e^) = 
0{e^). So the magnetic contribution to 5n|bdry is indeed of order and thus subleading when compared with the 
rhs of (4.29), that is of order e. 
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o)R0 




Figure 1: Plot of the dimensionless dispersion relation (jj{k) for the Rayleigh-Plateau instability in a 
static uniform tube for several values of the axial magnetic field B(o) (in the plot, B). The instability 
strength and decreases as the magnetic field grows. We use (4.32) and the numerical data correspond to 

tension potential energy can decrease (to leading order in the perturbation), and thus the condition 
(4.33) for the instability is independent of the introduction of a background magnetic field With 
or without a non-dynamical magnetic field, the Rayleigh-Plateau instability is a long wavelength 
instability that afflicts the plasma tube when its length is bigger than its transverse radius. 

The presence of the non-dynamical magnetic field does however reduce the strength of the 
instability. Indeed, in the plot of the dispersion relation (4.32), see Fig. 1, we find that an increment 
in the magnetic field ^8(0) weakens the instability. For high values of i3(o) the instability strength 
becomes considerably weaker. Strictly speaking we cannot however say that it disappears for a 
critical value of ;B(o) since the threshold mode (4.33) is always non- vanishing independently of B{Qy 

A word of caution is in order. Prom the equation of state (4.16) and the equilibrium condi- 
tion (4.17) one gets for d = 4: pj-Q) + Pj-g) = 5ah{i',b) — 2;B(o)A^(o) (for constant a). Note that 
unfortunately wc do not know the function /i(z^, b) in five dimensions and this prevents an accurate 
computation. Depending on its form, the condition p(o) + ^(0) > might require a maximum 
allowed .B(o)- To plot Fig. 1 we took 5a/i(z^, b) = 1, i.e., p(o) +^(o) = 5 — 2;S(o)-M(o)) a-^d this indeed 
implies a critical -6(0). However, we should have in mind that h is not a constant but h = h{i', b). 
There might thus not exist a maximum B(o) ■ 

4.5 Unstable plasma tubes and their gravitational duals 

The magnetized plasma tubes we have been discussing are dual to magnetic black strings in SS 

compactified AdSe- The Rayleigh-Plateau (RP) instablity of the plasma tube is dual to the Gregory- 
Laflamme (GL) instability of the black strings. Unfortunately the magnetic black string solutions 
in SS compactified AdSe have not been constructed yet, and thus we cannot compare our plasma 
results with their gravitational duals. Therefore, our study provides solid predictions for the features 

'^^Later, in the end of next subsection, we will argue tliat we do no expect this to be true when the approximation 
5x ~ is relaxed. 
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and stability of SS AdSg magnetic black strings. 

There is however growing evidence that justifies the comparison of our plasma results with 
the properties of magnetically charged black strings in an asymptotically flat background. Let us 
pause here to review the evidence that seems to permit this extrapolation. In [2, 33] the phase 
diagram of axisymmetric rigidly rotating plasma configurations was found. The set of solutions is 
given by plasma balls, pinched balls and plasma rings. Quite remarkably, the entropy vs angular 
momentum (for fixed energy) phase diagram displays remarkable similarities with the phase diagram 
of asymptotically flat [40, 41] and AdS [42] black hole solutions. 

More closely connected with our present study, [6] analyzed the RP instability of a neutral 
plasma tube and the associated phase diagram of non-uniform plasma tubes and localized plasma 
balls. The plasma instability has properties remarkably similar to the GL instability of a neutral, 
asymptotically fiat black string [7] (sec [43,44] for reviews). A brief account of these similarities 
includes [6]: (i) the dispersion relation of the two instabilities are quite similar and, not less 
important, both evolve similarly as we add extra dimensions (compare Fig. 5 of [6] with Fig. 
3 of [45]); ii) for static configurations, only axisymmetric modes are unstable and there is even 
a quantitative match for the threshold wavenumber as the number of dimensions gets large; for 
rotating solutions, non-axisymmetric modes can become unstable in both instabilities; iii) there are 
similar critical dimensions for the entropically favored solutions; iv) the phase diagram containing 
the solutions living in a compact dimension, namely (non-)uniform extended configurations and 
localized configurations is astonishingly similar (compare the S{J) diagram of Fig. 1 of [6] with 
Fig. 3 of [46]); the known initial time evolution of the GL instability of an asymptotically flat black 
string is similar to the dynamical evolution of the RP instability (compare Fig. 4 of [47] with Fig. 
1 of [48]). General arguments suggest that in the limit of large number of dimensions the fluid 
description of asymptotically flat black holes should become more and more accurate. 

The just mentioned facts motivate us to compare our results for the RP instability of a magne- 
tized plasma tube with the know studies [45, 49, 50] of the GL instability of magnetically charged 
asymptotically flat black strings/branes. In [49,50] it is found that a magnetic charge weakens the 
strength of the GL instability. In view of the above evidence, this is in agreement with our plasma 
result: a magnetic fleld also decreases the strength of the RP instability of the plasma tube. 

Besides showing that the GL instability gets weaker as magnetic charge is added to the system, 
[45,49,50] have also been able to show that it actually disappears in the extremal case, when the 
charge has its maximum allowed value (for some systems it is absent even before the extreme state 
is reached [50]). This property leads to an interesting new relation. Indeed, a dynamical instability 
like GL is expected to be accompanied by a local thermodynamic instability. This motivated 
Gubscr and Mitra [51] to propose the conjecture that a black brane/string with a non-compact 
translational symmetry is classically stable if and only if it is locally thermodynamically stable. A 
detailed analysis of this conjecture was further done in [52]. In these studies it was indeed found 
that there is a change in the sign of the specific heat of the black brane precisely for the critical 
charge where the GL instability disappears. It is thus interesting to discuss a natural version of 
the Gubser-Mitra conjecture in the uniform plasma tube setup. Unfortunately, we are not able to 
address this interesting issue. The reason is two-folded. First, to compute the specific heat of the 
uniform plasma tube we would need to have the relation between the entropy and the energy of 
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the system. To get the entropy we would need to integrate the entropy density given by the first 
relation in (4.15). This requires knowing the precise form of the thermodynamic function h{i',b), 
defined in (4.9), in five dimensions. Unfortunately, this is currently unknown. Second, although 
the RP instability gets weaker as the magnetic field increases we do not find a finite critical value 
above which it is not active (see however the discussion at the end of subsection 4.4). Wc do find 
that the instability strength reduces substantially as B grows not too large. Although in practice 
we can say that it disappears for large magnetic field, strictly speaking we are not allowed to say 
that it is gone since the threshold wavenumber keeps finite (see Fig. 1). The RP modes become sort 
of marginally stable in the sense that unstable modes are present but with very small strength. So, 
at most we would be able (if we knew h{iy, b)) to check that the specific heat also never becomes 
positive for the uniform plasma tube as the magnetic field grows. In the uncharged case, /i(i^, b) is 
a constant, and wc found in [6] that the heat capacity is indeed always negative. Adding a small 
magnetic field should not change this, but we cannot explore the situation for larger fields. 

The absence of a critical magnetic field in the RP instability deserves however a much more 
careful look. As described above this occurs because the threshold wavenumber, where the RP is 
marginally stable, persists invariant under changes of the magnetic field. This is an odd feature^^ 
and it is not hard to convince ourselves that this is a consequence of neglecting the susceptibility 
perturbations (see discussion associated with (4.19) and footnote 17). Indeed, from footnote 17 
we see that when 6x "j^ there is a contribution proportional to the square of the magnetic field 
strength in the Ihs of (4.29). This contribution persists in the following equations and the threshold 
wavenumber (4.33) gets then reduced when the magnetic field increases. Unfortunately, when we 
consider Sx 7^ 0, the perturbed system is left undetermined even after using the Bianchi identities 
and the conservation of the polarization current (in footnote 2). We would need to know some 
microscopic information to express dx in terms of perturbations of other thermodynamic quantities 
- a relation that we do not have. For this reason, and because Sx is nevertheless expected to be 
small (as justified previously), we took (5% ~ in our analysis. This approximation nevertheless 
allowed to find that a magnetic field decreases the RP instability strength, and we have a good 
control on the effects that a non-vanishing susceptibility perturbation introduces in the system: the 
instability strength still decreases when B grows and, on the top of this, the critical wavenumber 
also decreases. 

Usually an unstable mode signals a bifurcation to a new branch of configurations in the phase 
diagram of solutions. Without magnetic fields this is indeed the case and the RP unstable mode of 
the uniform plasma tube leads to a branch of non-uniform plasma tubes that joins a third branch 
of plasma balls localized on the compact tube direction [6,54]. Although not studied here, again 
because we do not know the function h{i',b), we naturally expect that similar new branches of 
solutions, representing magnetic non- uniform tubes and magnetic localized plasma balls, exist also 
in the phase diagram of magnetized plasma configurations. 

To end this section, we address the effects that viscosity introduces in the Rayleigh-Platcau 
instability. These were already discussed in more detail in [6], so wc will be brief in our comments. 

^^Indeed, if we consider the standard RP instability on a charged fluid tube coupled to a dynamical magnetic field, 
we find that both the instability strength and the critical wavenumber decrease when the magnetic field increases [53] . 
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Generically, at least classicaly, the viscosity increases the wavelength of the most unstable mode, 
and weakens the strength of the instability [53]. It has a subleading effect on the activation of 
the instability (however, the instability can get considerably weaker if the viscosity is very high), 
but plays an important role at later stages in the time evolution of the instability, namely in the 
pinch-off phase of the plasma tube [48]. Typically, the lower the viscosity is, the higher is the 
number of plasma balls formed. Your fluid is non-conformal so we have bulk and shear viscosity 
and thermal dissipation. Dissipation increases considerably the technical challenge of solving (even 
numerically) the MHD perturbed equations: we get a coupled system of differential equations with 
a fourth order derivative term. 



4.6 Regime of validity 

To conclude this section, the regime of validity of our results should be kept in mind. 

First of all, the fluid description of the deconfined plasma must be accurate, and hence the 
thermodynamic quantities must vary slowly over the mean free path £^{p of the constituent particles, 
which is of the order of the mass gap of the theory, or equivalently of the order of the deconfinement 
temperature Tg ~ ^. That is, all length scales A in the fluid must be 

A » T"^ ~ — . (4.34) 
Po 

This condition imposes restrictions on the validity of the Rayleigh-Plateau instability analysis: the 
RP unstable frequencies and wavenumbers must satisfy 

{LoRo,kRo}:^ (4.35) 

PO-fM) 

Since the most unstable mode dominates the instability we must guarantee that this condition is 
verified in the vicinity of the maximum of the dispersion relation. We find that for ^^^^ ^ 10^'^ 
this relation is satisfied, and things get better as becomes smaller. In particular this is true 
in the dispersion relation plot of Fig. 1. 

Second, the interface between the confined and deconfined phases, i.e., the fluid surface, has a 
finite thickness of the order 1/Tc [32], and therefore the delta-like surface approximation we used 
is valid provided that the curvature of the surface is small with respect to the scale l/Tc- So, for 
plasma tubes the analysis is good when the boundary radius Rg is everywhere large when compared 

with 

(4.36) 

which is perfectly compatible with (4.35). 

Finally, we neglected the dependence of the surface tension on the temperature and other 
thermodynamic quantities of the fluid. For consistency we must then demand that on the boundary 
between the confined and deconfined phases the temperature of the plasma must remain everywhere 
close to the critical temperature Tc, i.e., T/Tc ~ 1 at the tube boundary. 

As discussed already in great detail in sTibsection 4.5, we also neglected the dependence of 
the magnetic susceptibility on the temperature and on other thermodynamic quantities. In the 
previous section we already identified the effects of taking this assumption and of relaxing it. 
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Once these conditions are fulfilled, the plasma approximation holds and can be trusted to 
study the properties of black holes. The moral of this analysis is that there is a broad window of 
parameters for which the analysis of the Rayleigh-Plateau instability is well within the required 
validity regime. 

Acknowledgements 

We thank Roberto Emparan for useful discussions. MC and OJCD are grateful to the organiz- 
ers of the workshop "Quantum Black Holes, Braneworlds and Holography", Valencia, May 2008; 
OJCD and DK thank CERN for hospitality during the programme "Black Holes: A Landscape 
of Theoretical Physics Problems", August-October 2008; OJCD thanks the Niels Bohr Institute 
for hospitality and the organizers of the workshop "Mathematical Aspects of General Relativity" , 
Copenhagen, April 2008, where part of this work was done. MC was supported in part by DURSI 
2005 SGR 00082, MEG FPA 2004-04582-G02 and FPA-2007-66665-G02, the European Community 
FP6 program MRTN-GT-2004-005104, the FWO - Vlaanderen, project G.0235.05 and in part by the 
Federal Office for Scientic, Technical and Cultural Affairs through the Interuniversity Attraction 
Poles Programme Belgian Science Policy P6/11-P. OJCD acknowledges financial support provided 
by the European Community through the Intra-European Marie Curie contract MEIF-CT-2006- 
038924. This work was partially funded by FCT through project PTDC/FIS/64175/2006. DK was 
supported in part by INFN, MIUR-PRIN contract 20075ATT78, and by the European Community 
FP6 program MRTN-CT-2004-005104. 

References 

[1] O. Aharony, S. S. Gubser, J. M. Maldaccna, H. Ooguri and Y. Oz, "Large N field theories, 
string theory and gravity," Phys. Rept. 323 (2000) 183 [arXiv:hep-th/9905111]. 

[2] S. Lahiri and S. Minwalla, "Plasmarings as dual black rings," arXiv:0705.3404 [hep-th]. 

[3] J. Evslin and C. Krishnan, "Metastable black saturns," JHEP 0809 (2008) 003 
[arXiv:0804.4575 [hep-th]]. 

[4] S. Bhattacharyya, V. E. Hubeny, S. Minwalla and M. Rangamani, "Nonlinear fluid dynamics 
from gravity," JHEP 0802 (2008) 045 [arXiv:0712.2456 [hep-th]]. 

[5] M. Van Raamsdonk, "Black hole dynamics from atmospheric science," JHEP 0805 (2008) 106 
[arXiv:0802.3224 [hep-th]]. 

[6] M. M. Caldarelli, O. J. C. Dias, R. Emparan and D. Klemm, "Black holes as lumps of fluid," 
arXiv:0811.2381 [hep-th]. 

[7] R. Gregory and R. Laflamme, "Black strings and p-branes are unstable," Phys. Rev. Lett. 70 
(1993) 2837 [arXiv:hep-th/9301052]. 



34 



[8] S. A. Hartnoll, C. P. Herzog and G. T. Horowitz, "Building an AdS/CFT superconductor," 
arXiv:0803.3295 [hep-th]. 

[9] S. S. Gubser, "Breaking an Abelian gauge symmetry near a black hole horizon," 
arXiv:0801.2977 [hcp-th]. 

[10] S. A. Hartnoll and P. Kovtun, "Hall conductivity from dyonic black holes," Phys. Rev. D 76 
(2007) 066001 [arXiv:0704.1160 [hep-th]]. 

[11] S. A. Hartnoll, P. K. Kovtun, M. Miillcr and S. Sachdev, "Theory of the Nernst effect near 
quantum phase transitions in condensed matter, and in dyonic black holes," Phys. Rev. B 76 
(2007) 144502 [arXiv:0706.3215 [cond-mat.str-el]]. 

[12] G. Compere and D. Marolf, "Setting the boundary free in AdS/CFT," Class. Quant. Grav. 
25 (2008) 195014 [arXiv: 0805. 1902 [hep-th]]. 

[13] S. Bhattacharyya, S. Lahiri, R. Loganayagam and S. Minwalla, "Large rotating AdS black 
holes from fluid mechanics," arXiv:0708.1770 [hep-th]. 

[14] J. Hansen and P. Kraus, "Nonlinear magnetohydrodynamics from gravity," arXiv:081 1.3468 
[hep-th]. 

[15] N. Andersson and G. L. Comer, "Relativistic fluid dynamics: Physics for many different 
scales," Living Rev. Rel. 10 (2005) 1 [arXiv:gr-qc/0605010]. 

[16] A. Ishibashi and R. M. Wald, "Dynamics in non-globally hyperbolic static spacetimes. Ill: 
anti-de Sitter spacetime," Class. Quant. Grav. 21 (2004) 2981 [arXiv:hep-th/0402184]. 

[17] D. Marolf and S. F. Ross, "Boundary conditions and new dualities: Vector fields in AdS/CFT," 
JHEP 0611 (2006) 085 [arXiv:hep-th/0606113]. 

[18] V. Balasubramanian and P. Kraus, "A stress tensor for anti-de Sitter gravity," Commun. Math. 
Phys. 208 (1999) 413 [arXiv:hep-th/9902121]. 

[19] M. Bianchi, D. Z. Freedman and K. Skenderis, "Holographic Renormalization," Nucl. Phys. B 
631 (2002) 159 [arXiv:hep-th/0112119]. 

[20] A. M. Awad and C. V. Johnson, "Holographic stress tensors for Kerr-AdS black holes," Phys. 
Rev. D 61 (2000) 084025 [arXiv:hep-th/9910040]. 

[21] D. T. Son and A. O. Starinets, "Viscosity, Black Holes, and Quantum Field Theory," Ann. 
Rev. Nucl. Part. Sci. 57 (2007) 95 [arXiv:0704.0240 [hep-th]]. 

[22] M. M. Caldarelli, G. Cognola and D. Klemm, "Thermodynamics of Kerr-Newman-AdS black 
holes and conformal field theories," Class. Quant. Grav. 17 (2000) 399 [arXiv:hep-th/9908022]. 

[23] S. W. Hawking, C. J. Hunter and M. Taylor, "Rotation and the AdS/CFT correspondence," 
Phys. Rev. D 59 (1999) 064005 [arXiv:hep-th/9811056]. 



35 



[24] S. Bhattacharyya, R. Loganayagam, S. Minwalla, S. Nampuri, S. P. Trivedi and S. R. Wadia, 
"Forced Fluid Dynamics from Gravity," arXiv:0806.0006 [hep-th]. 

[25] A. H. Chamseddine and W. A. Sabra, "Magnetic strings in five-dimensional gauged supergrav- 
ity theories," Phys. Lett. B 477 (2000) 329 [arXiv:hcp-th/9911195]. 

[26] D. Klemm and W. A. Sabra, "Supcrsymmctry of black strings in D = 5 gauged supergravities," 
Phys. Rev. D 62 (2000) 024003 [arXiv:hep-th/0001131]. 

[27] K. Copscy and G. T. Horowitz, "Gravity dual of gauge theory on S'^ x x M," JHEP 0606 
(2006) 021 [arXiv:hep-th/0602003]. 

[28] R. B. Mann, E. Radu and C. Stelea, "Black string solutions with negative cosmological con- 
stant," JHEP 0609 (2006) 073 [arXiv:hep-th/0604205]. 

[29] Y. Brihaye, E. Radu and C. Stelea, "Black strings with negative cosmological constant: 
Inclusion of electric charge and rotation," Class. Quant. Grav. 24 (2007) 4839 [arXiv:hep- 
th/0703046]. 

[30] A. Bernamonti, M. M. Caldarelli, D. Klemm, R. Olea, C. Sieg and E. Zorzan, "Black strings 
in AdSs," JHEP 0801 (2008) 061 [arXiv:0708.2402 [hep-th]]. 

[31] T. Delsate, "New stable phase of non uniform black strings in AdS^," arXiv:0808.2752 [hep-th]. 

[32] O. Aharony, S. Minwalla and T. Wiseman, "Plasma-balls in large N gauge theories and local- 
ized black holes," Class. Quant. Grav. 23 (2006) 2171 [arXiv:hep-th/0507219]. 

[33] S. Bhardwaj and J. Bhattacharya, "Thermodynamics of Plasmaballs and Plasmarings in 3-1-1 
Dimensions," arXiv:0806.1897 [hep-th]. 

[34] V. Cardoso and O. J. C. Dias, "Gregory-Lafiamme and Rayleigh-Plateau instabilities," Phys. 
Rev. Lett. 96 (2006) 181601 [arXiv:hep-th/0602017]. 

[35] V. Cardoso and L. Gualtieri, "Equilibrium configurations of fluids and their stability in higher 
dimensions," Class. Quant. Grav. 23 (2006) 7151 [arXiv:hep-th/06 10004]. 

[36] V. Cardoso, O. J. C. Dias and L. Gualtieri, "The return of the membrane paradigm? Black 
holes and strings in the water tap," Int. J. Mod. Phys. D 17 (2008) 505 [arXiv:0705.2777 
[hep-th]]. 

[37] O. J. C. Dias, T. Harmark, R. C. Myers and N. A. Obers, "Multi-black hole configurations on 
the cylinder," Phys. Rev. D 76, 104025 (2007) [arXiv: 0706. 3645 [hep-th]]. 

[38] U. Miyamoto and K. i. Maeda, "Liquid bridges and black strings in higher dimensions," Phys. 
Lett. B 664 (2008) 103 [arXiv:0803.3037 [hep-th]]. 

[39] U. Miyamoto, "Curvature driven diffusion, Rayleigh-Plateau, and Gregory-Laflamme," Phys. 
Rev. D 78 (2008) 026001 [arXiv:0804.1723 [hep-th]]. 



36 



[40] R. Emparan and H. S. Reall, "A rotating black ring in five dimensions," Phys. Rev. Lett. 88 
(2002) 101101 [arXiv:hep-th/01 10260]. 

[41] R. Emparan, T. Harmark, V. Niarchos, N. A. Obcrs and M. J. Rodriguez, "The Phase 
Structure of Higher-Dimensional Black Rings and Black Holes," JHEP 0710 (2007) 110 
[arXiv:0708.2181 [hep-th]]. 

[42] M. M. Caldarclli, R. Emparan and M. J. Rodriguez, "Black Rings in (Anti)-deSitter space," 
arXiv:0806.1954 [hep-th]. 

[43] B. Kol, "The phase transition between caged black holes and black strings: A review," Phys. 
Rept. 422 (2006) 119 [arXiv:hep-th/0411240]. 

[44] T. Harmark, V. Niarchos and N. A. Obcrs, "Instabilities of black strings and branes," Class. 
Quant. Grav. 24 (2007) Rl [arXiv:hep-th/0701022]. 

[45] R. Gregory and R. Laflamme, "The Instability of charged black strings and p-branes," Nucl. 
Phys. B 428 (1994) 399 [arXiv:hep-th/9404071]. 

[46] T. Harmark and N. A. Obers, "Phases of Kaluza-Klein black holes: A brief review," arXiv:hep- 
th/0503020. 

[47] M. W. Choptuik, L. Lehner, I. Olabarrieta, R. Petryk, F. Pretorius and H. Villegas, "To- 
wards the final fate of an unstable black string," Phys. Rev. D 68 (2003) 044001 [arXiv:gr- 
qc/0304085]. 

[48] M. Tjahjadi, H. A. Stone and J. M. Ottino, "Satellite and subsatellite formation in cappilary 
breakup," J. Fluid Mech. 243 (1992) 297. 

[49] R. Gregory and R. Laflamme, "Evidence For Stability Of Extremal Black P-Branes," Phys. 
Rev. D 51 (1995) 305 [arXiv:hep-th/9410050]. 

[50] T. Hirayama, G. w. Kang and Y. o. Lee, "Classical stability of charged black branes and the 
Gubser-Mitra conjecture," Phys. Rev. D 67 (2003) 024007 [arXiv:hep-th/0209181]. 

[51] S. S. Gubser and I. Mitra, "Instability of charged black holes in anti-de Sitter space," arXiv:hep- 
th/0009126. 

[52] H. S. Reall, "Classical and thermodynamic stability of black branes," Phys. Rev. D 64 (2001) 
044005 [arXiv:hep-th/0104071]. 

[53] S. Chandrasekhar, "Hydrodynamic and hydromagnetic stability" (Dover Publications, New 
York, 1981). 

[54] K. i. Maeda and U. Miyamoto, "Black hole-black string phase transitions from hydrodynam- 
ics," arXiv:081 1.2305 [hep-th]. 



37 



